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QJ , Bogolubov's chain of equations method for temperature Wightman functions is suggested for 

[T i . investigation of relativistic phase transition. The chain equations for the Wightman func- 

tions forming momentum-energy tensor are obtained. It is clarified that structure of the 
chain equations determines the basis approximation (the Hartree - Fock approximation) and 
corrections calculation algorithm. The basis approximation is investigated in details: renor- 
malized equations for effective masses, order parameter and generating functional which 
reproduce those equations are obtained. Being considered on the solution of the gap equa- 
tions for the effective masses the generating functional turns to nonequilibrium functional 
of free energy density, which allows to obtain phases stability conditions. Thermodynamic 
observables like heat capacity and sonic speed are calculated. The correction to the Hartree- 
Fock approximations is ascertained to be small for all temperatures excluding vicinity of the 
phases equilibrium point. 
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1 Introduction. 
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An idea that phase transition may take place in a quantum field system was first proposed by 



Kirznits|T[. First quantitative calculations were made by Kirznits and Linde 0, and later by 
Dolan and Jackiw M, Weinberg M. Theory of relativistic phase transition (RPT) deals with 
nonequilibrium Landau functional (which is usually called "effective potential"). When quasi- 
^ particle approximation is used for calculation of effective potential from microscopic operational 
Lagrangian, a problem of quasi-particles mass spectra arises. Calculations made in Matsubara 
diagram technique showed, that phase transitions can't be described in iterative perturbation 
theory approach — tachion pole in mass spectra of scalar particles and consequently imaginary 
terms in Matsubara "effective potential" appears far from the point of phase transition. En- 
hanced resuming schemes successfully solved the problem of tachion pole in mass spectra but 
not the problem of effective potential calculation. 

As for calculation of nonequilibrium functional in Matsubara approach, with temperature 
Green functions obtained in the framework of resuming schemes, the extremum condition of 
such functional leads to equation for order parameter which is different from that obtained 
directly from equations of motion. Illusion of solving of the problem mentioned is achieved 
only by manipulation with high-temperature approximation of Bose-Einstein integrals. 

Cornwall, Jackiw and Tomboulis suggested an "effective functional" of different kind |J. In 
addition to equation for order parameter the " C JT-functional" reproduces equations for full 
Green functions. At present O(N)— model is almost completely investigated in this approach 
in Hartree-Fock approximation. Nevertheless problems of renormalization of the functional, 
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analytical analysis of the phases stability conditions, calculation of observables and estimation 
of Hartree-Fock approximation accuracy need to be discussed in details. 

These problems (at least for O(N)— model) can be solved in alternative approach, based on 
Bogolubov's chain of equations for temperature Wightman functions (WF). 

In Wightman approach the equations of motion and observables of a quantum-field system 
are written in terms of Wightman functions. || . In this paper equations for all the functions 
forming energy-momentum tensor are presented. Field operators smoothing, used in axiomatic 
Wightman formalism, replaced with scheme of dimensional regularization. Self-consistency of 
that redefinition method is proved. 

It is found out that effective masses of quasi-particles, which are formed by particles inter- 
actions with each other and vacuum, automatically arise in the Bogolubov's chain equations. 
Thus, it is natural to use Hartree-Fock approximation as a basis approximation of reduction. 
Corrections calculation algorithm is also determined by the Bogolubov's chain mathematical 
structure. 

Hartree-Fock approximation was investigated in details. A functional (generating func- 
tional) which reproduce gap equations for effective masses and equation for order parameter is 
obtained. It is shown, that the generating functional can also be calculated from momentum- 
energy tensor components, which proves self-consistency of thermodynamic equations of state 
obtained from quantum equations of motion with those obtained from generating functional. 
The generating functional, being considered on solutions of the gap equations for effective 
masses turns into functional of nonequilibrium free energy, which is used for analytical analysis 
of the phases stability condition. Functional of equilibrium free energy, which is obtained by 
substitution of the order parameter on temperature dependence into nonequilibrium one, is 
used for calculation of thermodynamical observables: heat capacity and sonic speed. 

An important problem of the Hartree-Fock approximation validity was clarified by calcula- 
tion of correction of the first order of vanishing to this approximation. 

A Lagrangian of the O(N)— model and equations of motions are written out in chapter ^|. A 
definition of Wightman functions, symmetrized with regard to permutation of it's arguments, 
chain's equations for Wightman functions up to fourth rank and chain reduction algorithm 
are given in chapter |3|. Equations for effective masses and order parameter, functionals and 
observables in Hartree-Fock approximation, and results of numerical calculations are given in 
chapter |j. Finally, chapter [5] is devoted to calculation of the first order of vanishing corrections 
to effective masses, order parameter and observables. 



2 Lagrangian and equations of motion. 

Operator equations of motions for quantum fields ore obtained from Lagrangian of system of 
N scalar fields <p a 

C = \ (d,<Pa ■ d^ a + /i 2 a • <j> a ) - A(0 O . <f> a ) 2 , (1) 
by standard variational procedure: 

4A 

<9 M <9 M a - //(pa + — (060fe)0 a = 0. 

For /i 2 < 0, the theory is invariant under O(N) transformations of the quantum fields. For 
fi 2 > and low temperatures spontaneous symmetry breaking described by nonzero vacuum 
expectation value of (f)^, 
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breaks down symmetry group to 0(N — 1) and leads to iV — 1 Goldstone bosons. 

(f) N = (p + V, 

<t>a = Xa, a = 1 . . . N - 1 



Equations for quantum fields Xa and <p look like: 

4A 

d^Xa - fJ?Xa + -J^( V Xa + 2v<pXa + v\a + XbXbXa) = 

Momentum-energy tensor is obtained from Lagrangian (|l|) by metric variation procedure. 
For homogeneous, isotropic and stationary system, momentum-energy tensor turns to: 



(2) 



(3) 



/j, 2 v 2 



Xv 4 



+ 



A 



N N 



(T;> = (d^if) + (d,Xad V Xa)- 

(<P 4 ) + (XaXaXbXb) + 2(lfi 2 XaXa) + 2v(ip 3 ) + 2v((fXaXa) 



(4) 



3 Bogolubov's chain of equations. 

Definition of symmetrical Wightman functions. Dynamics of a system, which is deter- 
mined by Lagrange equations of motions, can be described in terms of chain of equations for 
Wightman functions (WF) — so called Bogolubov's chain. 

Full n— point WF (WF of rank n) is an over state expectation value of n quantum field 
operators, taken, generally speaking, at different space-time points. 

{f{x 1 )^P{x 2 )-f{xi)-X{x m )-X{x n ) )■ 

It can be represented as a sum of various WF of lower rank (the sum of ranks in each product 
is equal to n) and so called correlative Wightman function of the same rank. This representa- 
tion, essentially, is a definition of correlative WF. For homogeneous, isotropic space-time it is 
convenient to use WF symmetrized with regard to permutation of it's arguments. Algorithm 
of it's definition, which takes into account residual 0(N — 1) symmetry of the theory under 
consideration, is illustrated for two-point G( XlX2 ), D( XlX2 ) and four-point G 22 ( XlX2 \ X3X4 ) func- 
tions. Here left bottom index is the number of (p fields and the right one is the number of Xa 
fields which form the corresponding Wightman function. 

1. An expectation value of a sum of every possible permutation of quantum-field operators 
is taken (symmetrization operation): 

2G(x 1 x 2 ) = ( i P(xi)^(x2))sym = (f(xi)f(x 2 )) + ( f(x 2 )¥(xi ) ) , 
2(N - l)D {xiX2 ) = (Xa(x 1 )Xa(x 2 )) sym = (Xa( Xl )Xa(x 2 )) + ( Xa(x 2 )Xa(x 1 ) ) , 
A(N -l)G 22{xiX2 \ X3Xi) = (f(x 1 )f(x 2 )Xa(x 3 )Xa(x 4 )) sym = ((P(x 1 )^P{x 2 )Xa(x s )Xa(x 4 ))+ 
+ (V(x 1 )V(x 2 )Xa(x 4 )Xa(x 3 )) + ( ¥(x 2 )f{ Xl )Xa(x 3 )Xa{x 4 ) ) + ( ( P(x 2 ) ( f(x 1 )Xa(x 4 )Xa(x 3 ))- 

Since WF should be 0(N — 1) invariant, quantum-field operators Xa come in pairs with the 
same indexes (Einstein sum rule is used). 

In that way defined functions are invariant of permutation of arguments of (p operators and 
of Xa operators with the same indexes. 
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2. Four-point WF G 22 ( XlX2 \ X3X4 ) can be represented as sum of various products of (non- 
symmetrized) WF of the second rank and corresponding correlative function: 

4 (AT — l)G^22(;ri X2 \xa X 4 ) ~ Xa(x 3 )Xa(x 4 Xa(x 4 ) Xa(x3 ) ) + 

( ( P(x2) ( P(x 1 ))(Xa(x 3 )Xa(x 4 ) ) + ( ^(x 2 ) ( P{x 1 ) ) ( Xa(x 4 )Xa(x 3 ) ) + 4(N - l)C 22 ( Xl X2 \ X3 X4 ) 

3. The sum terms can be regrouped, which picks out symmetrized WF: 



G 



22(X!X2 \x 3 X 4 ) 



G( Xl x2)D(x 3 x 4 ,) ~\~ C 22 ( XlX2 \ X3XA y 



It should be noted that permutation of arguments separated by vertical line leave a Wightman 
function, in particular G 22 ( Xl X2 \ X3X4 ) } invariant: G 22 ( XlX2 \ X3X4 ) = G 22 ( X2Xl \ X3Xi ) = G 22 ( XlX2 \ X4X3 } = 

^22(x 2 x 1 \x 4 x 3 )- 

It is convenient to divide WF which depend on odd number of arguments by expectation 
value v. 

Residual global and discrete (x a ~Xa) symmetries of the theory restrict number of 
nontrivial WF. In particular, there exists only two nontrivial three-point 



G 03 (xi X2 x 3 ) 



1 



3! v 



sym 



r _J v 

^21 Oi X2 \x 3 ) — 2\ v N — I X- a ( x i)^ a ( x i) { P( x i ) / 



sym 



and three nontrivial four-point correlative WF: C 22 r XlX2 \ X3X4 ) 

1 



and 



G 04 (xi X2 x 3 x 4 ) ^| (f(xi) ( p(x2) l P(x 3 ) l P(x 4 )) sym G( Xl X2 )G( X3 Xi ) 

~G( Xl X3 )G[ X2 Xi ) — G( X1 X4 )G( X2 X3 ) 

/jy _ ^2 (Xa(x 1 )Xa(x 2 )Xb(x 3 )Xb(x 4 )) sym~ 



G 40 (xi X2 \x 



3x4) 



(N - 1)D( X1 X2 )D( X3 Xi ) + D( X1 Xz )D( X2 XA ) + D( X1 Xi )D {X2 X3 ) 



N — 1 

Correlative functions of higher rank are defined in the same way. Strictly speaking, written 
above Wightman functions are divergent and should be redefined by procedures of regularization 
and renormalization. Q In this paper the scheme of dimensional regularization is used. This 
technique lets retain exact thermodynamic relations between equations of state, functionals 
and observables after renormalization. 

In terms of Wightman functions momentum-energy tensor (Hf) looks like: 



(t;> = lim (d, {x) d» {xi) G {xxi) + (N — l)d, {x) dl Xl) D (xxi) ) - <y 



+ 



A 

+ N 



(xx) 



N 



Gq 4 (x X X x)~^ \P) 



+ (N — 1) 2 C 4[) („„) + 2(N — 1)C 22 („„) + 2v 2 C 0S ( xxx ) + 2v 2 (N — l)C 21 ( xx \ x ) 



1 In discussed formalism these operations replace quantum-field operators product smoothing, used in ax- 
iomatic quantum-field theory. 
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Derivation of the chain equations. First of all equations for WF forming momentum- 
energy tensor should be obtained. 

Equation of state for vacuum expectation value is obtained by averaging (0): 



4A 



-/i 2 + -jy (V + 3G {xx) + C 03[xxx) + (N - l)D (xx) + (N- l)C, 



21 (sX \x) 



0. (6) 



Averaging of ([|) gives 



2«< VXa ) + ( ¥ 2 Xa ) + ( XbXbXa ) = 0, 



which is identity owing to mentioned above symmetry Xa ~Xa- 

Averaging equations (|2|) and ([|) multiplied by one or several operators, one obtain equations 
of the Bogolubov's chain. After symmetrization, which is required in account of symmetrical 
definition of WF, commutators arise: 

A(f(xxi) = <f(x)<f(xi) - Pixi^ix), AX(xxi) = N _ 7 (Xa(x)Xa(x 1 ) ~ Xa(xi)Xa(x)) ■ 



Equations for two-point WF look like: 



d^G^ xxl ) + 



+ 



4A 

iV 



4A / 

-/i 2 + — (3v 2 + 3G (xx) + (N-l)D 



N 



(xx) 



G( XXl )+ 



{C 0i ( XXXXl ) + (N — l)C 22 ( X x\xxi)j + v 2 (3C 0S ( XXXl )+ (N — 1)C 21 ( XX \ X1 )j 



(7) 



d^D {xxi) + 
4A r 



4A / 

-/i 2 + — (v 2 + G (xx) + (N + l)D {xx 



N 



+ 



N 



(A^ l)C*4o(x a; |x xi ) G 2 2(xx\ \xx ) "I - 2l7 C 2 i (x xi \x ) 



D(xX! ) + 





Owing to homogeneity, isotropy and stationarity of the system under consideration G( X xi) an d 
D( XX1 ) depend only on module of two arguments difference \x — X\ |; consequently, when x = X\ 
their value are determined only by internal system parameters. Thus, emerge in equations (|TTD, 
(|T2]), all equations for WF and equations for commutators values 



4A 



m\ = - fi 2 + — (3v 2 + 3G (xx) + (A — 1)D { 



xx) 



4A 



m 2 2 = -/i 2 + — [v 2 + G [xx) + (N + l)D {xx) 



(9) 



(10) 



depend on the system state but not the 4-coordinates. Hence these values play role of param- 
eters when chain equations are solved. 

Taking into account the designations introduced, equations for two-point WF can be rewrit- 
ten as follows: 



+ 



4A r 

A 



d^d^G^xi) +ml G( XXl )+ 

(XXX X\ 

) + (A - l)C 22 ( xx \ XXl )^j + v 2 (3C 03 ( XXXl )+ (N - l)C 21 ( xx \ xi y) 



+ 



4A 

iV 



dpdPD^) +m\ D {xxi) + 

(A 1 )C 40 ( X x\xxi) C22 (x xi \x x ) 2f G 21 ( XX1 \ x ) 
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0. 



(11) 



(12) 



Equation for WF C 03 („ ll2 ) looks like: 

+ iV ^ (xa:i) ( 3 ^( IJ!I ») + (^V- l)C a i(**|a:a)) + G( XX2 ) (^C m ( XXXl ) + 

4" (-^ 1)^21 {XX \xi_ ) J "I - (-^ l)^22(ffi £ |£l £2 ) 4~ 3Cq4(z Mi 12 ) "I - Cos(£ £ £ £l £2 ) 4~ 

24A U 1 A 



(13) 



iV 



+ (iV 1)C 23 ( 3 ; 3 ; a; 2 ) 

Since C 03 ( XX1X2 ) is invariant under transformations x <-> Xi , x <-> X2 and Xi <-> X2 , solution 
this WF is a sum of solution of ([13]) and solutions obtained from that by permutations 
-> Xi and x <-> x? . 



for this WF 
x <-> Xi and x <-> x 2 . 

For C 21 (x Xl \ X2 ), which depends on ip and Xa quantum fields, one should derive two equations, 
which describe contributions of <p and Xa respectively. 
The first one is derived from (|3]), 

d pd^G 21 ( x Xl \ X2 ) ^21 ^1 \x2 ) Jr 

4A r 

"I - 2C 22 (x x\xi X2) 4~ ^23(x x\x X \ X2) 4~ ^03 ("12) "I - X2 )C 21 (x xi \x)~^~ ( 14) 

(a; x \x xi \x 2 ) (x x \x 2 ) 

and the second one from @, 



dpd X2 C 21 ( XXl \ X2 ) + 77i 1 Cjj^jpj,^ |s 2 )~l~ 



4A r 

-^r [3C 22(:r:ri 1^^) + 1^ I2a;2 ) + (iV - 1)C 41 ) + 

4" (-AT l)C4o(x£i \x2 X2 ) 4~ 2AZ)( I:L X2 )C 21 (a;a;2 |asa ) 4~ 2AG( rr2 . 2 )C 21 (x 1 £ 2 |z 2 ) 



(15) 



BA 

'iV 



D(xi X2 )D( XX2 ) 4" ^X(£2 X )^X{X2 X\ ) 



Since C 21 r xxl \ X2 \ is invariant under transformation x <-> Xi solution for this WF is a sum of 

and solution of ( 14]) after permutation x <-> Xi . The same assertions 



solutions of (|i~5|) and ( |14] , , 
are valid for all the other WF 



Equation for C 04 ( XXlX2X3 ) is derived from 



4A r 



+ N 



G( XXl ) 3(7, 



d ll d^C 0i [ xxlX2X , i ^ + 777.J C 04 ( a;a;i3:2I3 ) + 

) 4" (-/V" 1)C* 22 (££ |£ 2 £ 3 ) +G( XX2 ) 3C Q4 ( XX2:3 Xl ) + 



y 04(£££2 £3 



+(A^ l)C 22 (a; x \ Xl X3 ) +G( a . X 3) 3C 04 ( XXXl X2 ) + (iV X)C 22 [ XX \ Xl X2 ) 
G( XXl )C 03 [ XX2X3 ) + G( xx . 2 )C 03 ( XXl x 3 ) + G( XX3 )C Q3 ( XXl X2 ) 



+ 



(ass I S3)) + Co3(££i £3 )( 3C 3 ( XXX 2 ) 

+ 



(16) 



+ u Co3(££i£2)(3C 3( ) + 

+ (-/V — 1)C 21 (££|£ 2 )) + C„3(££2£3)(3Co3(£££ 1 ) + (N ~ l) C 2X (3. jj. | a!1 ) ) 

05<£££1 £2 £3 ) (-^ 1)C23(£ £ |£1 £2 £3 ) ^06 (£ £ £ £l £2 £3 ) 



+ 



3C, 
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+ {N- l)C 2i{x x \xxi x 2 x 3 ) 



A. 

N 



C 2^G{xx 1 )G(xx 2 )G{xx 3 ) + 



+^ t P(xx 1 )^ ( P(xx 2 )G( XX3 ) + Acp ( x X2 )A<p( XX3 X1 ) + A<p( xxl -)A(p( XX3 ) G( XX2 ) 

Since the chain equations are quite bulky, a simple criterion of it's propriety is very useful. 
For example WF C 04 ( XXlX2X3 ) is invariant of permutation of it's arguments, so that equation 
(0) should be invariant under transformations x\ «-> x 2 , x\ x 3 , x 2 <-> x 3 (but, due to 
differential operator d^d%, not the permutations involve x). Easy to check that the equation, 
written above, satisfy this criterion. 

Equation for WF C i0 ( XXlX2X3 ) is derived from (|]): 



4A 
+ N 



d[id x C i0 ( xxl \ X2X3 } + w 2 G 40 ( xxlX2X3 ^-\- 

^ (ni)((-^ ~l" 1)^40 (as xx 2 x 3 ) G 2 2(x2 X3 \xx ) 4" 2f G 21 ( X1 X3 \ x )) + D( XX2 } (C 40 ( x x \ XlX3 )~\~ 



~\~2C w ( XXl \ XX3 ) + ^ G 22 ( X1 X3 \xx) + jy _ 1 ^21 (xi X3 |x)) + D( XX3 ) (C 40 ( I3 ; Xl 

~^^2l(a; X3 \x)G 21 (xi X 2 \x) 1)C 21 (jjjj I x )C 21 ^3 |x ) 4" 2t> C 41 (xijij 13 Ix) - ! - 

8A(iV + 1) 



+ C 42 ( XXlX2X3 \ xx) + (N - l)C 6o( 

111111213) 



N • (N — 1) 



D( XX3 )D^ XX2 ")D( XX1 y 



A 



iV(JV- 1) 



(5iV - 7) 



^X(x X2 )^X(x X 3 )D(x X\ ) "I - A\(^ XX1 ^A"%(xX2 )D(xX3 )~^~ 



+Ax(xx 1 )Ax(xx 3 )D( x 



x 2 ) 



For WF C 22 ( X2X3 \ X xi) one should obtain two independent equations. The first one is derived 
from (0) by multiplying by Xa{x 2 )Xa{x 3 )<f(xi ) and Xa{x 3 )Xa{x 2 )<f{xi): 



4A r 

iV 

~\~Di x ,„ \C, 



9fj,d x C 22 (x 2 x3 \xx 1 ) + m i G 22 (x 2X3 \xx\ ) + 
G( x x\ ) (i-N ^)G 40 ^ X xx2X3) 3C 22 ( X2X3 |xx)J "I" 2 D( XX2 )G 22 ( XX3 \ XXl ) + 



(113 ) L/ 22(xx2 |a; xi ) 



~\- 2, V 3G( X Xl )C 21 (x 2 x 3 \x) "I - -^(x X3 )G 21 (x X2\xi) D( x X2 )C 2 1 (x X3 |xi ) 



+3f C 2 i(s2S3 |x)C()3(xxxi ) "I" v (N ^)G 21 ( X2 X3 \ x )C 21 (xx |xi ) "I" 2C 21 (a; X3 \ Xl )C 21 (x X2 \ x )~\~ 



~^^G 2 ,f XX2 \ x , )C. 



(XX 2 \xi ) L/ 2l(«3 |x) 



-\- 3v C 23 ( X2 X3 \ x x Xl ^ ~\~ C 24 ( X2 X3 \ x x x Xl } ~\~ i^N l)f C 41 (x x X2 X3 1 2; 1 ) 



+ (iV l)C 42 (a;xx2X3 |xxj ) 



8A 

TV 



D(xx 2 )D(xx 3 )G(xx 1 ) + —Ax(xx 2 )Ax( X x 3 )G( XXl ) 



and the second one from (Q) by transformation x —>■ x 2 and multiplying by xi x 3 ) l f{x ) l fi(xi 
and x^M^i )<p(x): 



4A 



(X2 X 3 ) 



d/j.dx 2 C 22 (x 2 X3 \xxi) ~l~ m 2 G 22 ( X2 X3 \ XXl )~\~ 

G QA (xx\ X2X2 ) "I" ("^ "I" l)C22(X2 x 3 \xx) "I" 2 G( XX2 )C 22 (x 2 X3 jjjjj X2 ) + 
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~\~G(xi x 2 )C 22 ( X2 X3 \xx 2 ) + 2f D(x 2 x 3 )C m (xx 1 x 2 ) + G(xi x 2 )C 21 (x 2 x 3 |as) + 



+G(xx2 )^2l(*2 ^3 \xi ) 



^^2i(a;2 i'3 \x 2 )C"o3{xxi x 2 ) ~\~ C 2 i(x 2 x 3 \xi )G 03 (xx 2 x 2 



+ a 



21 (x 2 X 3 \x)C 03 (xi X2 Xa) 



+ t; 2 (iV+l) 



^21(22 X2 |xi )^2l(^2 ^3 |xi ) + ^21(^2 212 \xi )C 2 l{x 2 X 3 \x) 



+ 



+(N-1)C 4 
_8A 



(X2 X2 X 2 X3 I 



x xi ) ^ 2 i{x 2 X3 |x xi X2 X2 ) "I - 2f C23(x2 x 3 \x x\ x 2 ) 



1 



D(x 2 x 3 )G( X x 2 )G(x lX2 ) + —A{p(x 2X )A(p(x 2X1 )D{x 2 x 3 ) 



Equations for commutators, which present almost in all chain equations, look like: 

d IM d£A(p {x x 1 ) + m\ A(p { x Xl ) = 0, 
d^Ax(x Xl ) +ml Ax(xxO = 0. 



(17) 
(18) 



In high symmetry phase, where vacuum expectation value is equal to zero, symmetry with 
regard to ip — > —ip transformation is restored, consequently all WF which depend on odd 
number of <p quantum-field operators are identically equal to zero. The chain equations for 
this phase are easily derived from those for low-symmetry phase. 

The chain reduction. Two approaches to the chain reduction are known: model-approximative 
and iterative. 

In model-approximative approach a few first chain equations are solved, and the higher 
rank Wightman functions are evaluated as a certain combination of lower-rank WF according 
to an algorithm specified. Use of this approach in the considered formalism is difficult due 
to necessity of WF redefinition: in order to renormalize divergent integrals one should know 
energy spectrum, which, in turn, can be obtained only after renormalization. 

In iterative approach corrections to a basis approximation are calculated. Choice of basis 
approximation is dictating by the chain mathematical structure. 

A chain equation is a generalized D'Alamber equation 

<9 M <9£C nm (x...x m ) + m 2 C nm (x...x m ) = f {x ■ ■ .x m ) 

with combinations of correlative functions and field commutators as sources, which solution 
is a sum of general solution of the corresponding homogeneous equation and partial solution 
concerned with the sources. 

Solution of homogeneous equation, considered in chapter ||, should be taken into account 
only for two-point WF G( Xl X2 ) and D( Xl X2 ), which form effective masses mi, m 2 . For correlative 
functions, which describe weak effects of many-particles interaction, in considered iterative 
approach one should take into account only solutions concerned with the sources. Since sources 
contain small constant of self-interaction, correlative WF C 03 ( XlX2X3 ), C 21 ( XlX2 \x 3 ) are of the 
first order of vanishing in compare with Gi yXlX2 ) and D/ xlX2 y 

Independently of choice of basis approximation the following algorithm of calculations can 
be formulated: 

1.1. Derive solutions of the chain equations (considered in a certain approximation) keeping 
mi, m 2 and v as parameters. As a result one obtain WF as functions of coordinates, effective 
masses, order parameter and temperature: 



(x X\ ...X-nX-rri) 



W(xxi . . .x n x m ,mi,m 2 ,v,T) 



8 



1.2. Calculate WF at coincident points. Owing to homogeneity, isotropy and stationarity 
of the system under consideration, the result depends on temperature, m 1 , m 2 and v, but not 
4-coordinates: 



XX ...XX 



) = W(mi, m 2 , v,T) 



1.3. Substitute Wightman functions, calculated at coincident points, in (0), (0) and @ and 
solve this system of nonlinear equations for v, mi, rri2- After substitution of solutions 



one obtain WF as functions of coordinates and temperature. 

1.4. By known dependence WF on temperature calculate thermodynamic observables using 



The chain reduction basis approximation corresponds to allowing only for two-point G( XXl ) 
and D( XXl ) Wightman functions in equations (|II|) and (^), equation of state for order param- 
eter (^D and momentum-energy tensor (|5|). 

When calculating corrections to the basis approximation, which assumes calculation of 
higher rank WF, already known two-point Wightman functions are used: 

11. 1. Substitute G ( XlX2 ) and D ( XlX2 ), calculated in basis approximation, in equations for 
correlative WF, in particular C 03 ( Xl X2 X3 ) and C 21 ( Xl X2 \ X3 ), and obtain partial solution, concerned 
with sources. 

11. 2. Substitute obtained solutions for correlative WF into the chain equations for Gr xlX2 ), 
D( XlX2 ), obtain their solutions and so on... 

11. 3. Calculate v, m 2 and thermodynamic observables on temperature dependence 
(according to points 1.2 — 1.4 of algorithm). 

It is easy to see, that corrections are of 0(iV _1 ) order of vanishing, thus at the limit N — > oo 
the basis approximation becomes exact solution of the chain equations. 

4 Hartree— Fock approximation. 

The Bogolubov's chain basis reduction approximation is called Hartree-Fock approximation 



The chain equations in Hartree— Fock approximation. According to algorithm, stated 
above, at first step equations for two-point functions, with ni\ and rri2 considered as parameters, 
are written out: 



These equations are similar to those for noninteracting fields, thus it is natural to represent (p 
and Xa as follows: 



mi = toi(T), m 2 = m2(T),v = v(T) 



momentum-energy tensor (|5|). 



(HFA)| 



d^G o{xxi) +m 



I G o(xxi) = 0,d^D o(xXl) + m 2 2 D o{xix) = 0. 




1 



2 =p 2 + mi, 




(19) 



2 Which is often mentioned as "mean field approximation", "Hartree approximation". 
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with usual commutators for ap,bp^,ap,bp, a s operators. This representation corresponds to 
description of the system as an ideal gas of particles with masses dependent on temperature. 
At the second step two-point WF are calculated in coincident points: 



AG 



o{xx) 



1 



1 



V 



V 



XD 



(xx) 



p Vz x p exp — 1 
1 1 



p ' - 2 p ~p Ve 2P exp ^ - 1 



5 



p Ve lP/ 



p Ve 2 p 



(20) 



As it was already mentioned, divergent terms were renormalized in the framework of dimen- 
sional regularization method. Essential part of the method is that a theory is initially formu- 
lated in non-integer dimension space and transmutation parameters present in Lagrangian as 
factor multiplying dimensionless constants of interaction. Indeed, in Lagrangian ([I]), written 
for space-time of dimension D=l-\-n = l-\-3 — 2e, number of dimension of self-interaction 
constant A is L n_3 ; dimensional transmutation parameter is defined as A n = A/ n_3 . Analytic 
continuation to integer dimension space e — >0,_D^l + 3is made only after subtraction of 
s^ 1 poles (this operation is equivalent to redefinition of inoculating constants of the model). 
After renormalization of divergent terms ( p0[ ) one obtain: 



p 2Ve p 



A 



m 

16^2 



In 



m 2 \ 
A 2 /' 



A 



4tt/ 



-2 



oC-\ 



where C = 0.5772157 is Euler constant. 

At third step renormalized WF are substituted to (|), (0) and (||). 

It should be noted, that independently of renormalization method used, equation of state 
for order parameter (0) in Hartree-Fock approximation is written as: 



Gap equations for effective masses 
follows: 



2 $A n 

mi ~N V 







(21) 



and (|10|) for further purposes should be rewritten as 



Ji(mi,T) + 



m 



16tt 2 



1 In ^ 



A 2 



N 



IX(N 



(N-l) 



Ji(m 2 ,T) + 



^2 

16tt 2 



In 



m 2 

A 2 



(N + l)m 2 — (N — l)m 2 + 2/i 2 
iV(JV-l) 



i\(N + 2) 



3m 2 — m 1 + 2pf 



:22i 



where following designations for integrals over Bose-Einstein distribution are introduced: 



J n (m,T) 



2vr 2 



p dp 



VP 2 + m ' 2 ov „ yV+m 2 



n 



0,1,2, 



Temperature integrals satisfy recurrent relations: 



dJ n (m,T) 
dm 



-(2n - l)mJ n _i(m,T), 



dJ n (m,T) 2n 



m 



dT 



J„(m,T) + (2n - 1) — J„_i(m,T) 



Gap equation for effective mass at high-symmetry phase is obtained from (22) by zero filling 
of order parameter. It is easy to see, that at high-symmetry phase mi = m 2 = m. 
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It follows from ( pi]) and (p2) that in Hartree-Fock approximation the following effects are taken 
into account: 

1) vacuum influence on quasi-particles properties (order parameter v in definition of effective 
masses) . 

2) whole particles influence on condensate and properties of each quasi-particle (temperature 
integrals in equations for condensate and effective masses). 

3) zero-point oscillations of quantum fields influence on condensate and properties of each 
quasi-particle (renormalized vacuum integrals in equations for v and effective masses). 

Generating functional. Lagrangian ([!]) is a "generating functional" for operator equations 
of state, cause it reproduce these equations on it's extremals. Since Bogolubov's chain is equiv- 
alent to set of operator equations, a "generating functional" which reproduce chain equations 
must exist. The same assertion is valid for the chain reduced in any approximation. 

In self-consistent approximation this functional should reproduce gap equations for effective 
masses and equation of state for order parameter, or their linear combination. 

By equations (^1|) and (|22|) generating functional, which satisfy requirements 

0, hzH =0, h- =0. (23) 



\ dmi J \ dm^ , 

\ 1 / mi,v \ " / mi,v \ / mi,m2 

is reconstructed up to an arbitrary constant. It's explicit form: 

F(T,m 1 ,m 2 ,v) = -- [J 2 (m 1 ,T) + (N - l)J 2 (m 2 ,T)] + U(m 1 ,m 2 ,v), 

U( mi ,m 2 , v) = ^ In ^ ^ v , - H — ui j - (24) 



N 



(N + l)m\ + 3(N 



8X(N + 2) 
+^{m\ + (N — l)m 2 2 ) 




16A(iV + 2)' 



Gibbs potential. The goal of the theory of phase transitions is a calculation of transition 
temperature, phase co-existence region, and thermodynamic observables. The last can be 
obtained from momentum-energy tensor components (g), which for homogeneous and isotropic 
system is diagonal (T") = diag(e, — p, — p, — p). Since state of the system under consideration 
is defined by temperature and volume, one should use Gibbs potential for thermodynamic 
description. Due to chemical potential equality to zero, potential of free energy coincides with 
fi-potential, which is easily calculated via momentum-energy tensor spatial components: 

Q = - lim Pn V n = lim X ll n . (25) 

As (|25|) is divergent, it should be renormalized. In (§) two-point WF G( xx ) and D( xx ) are 
expressed through effective masses and, consequently, are not divergent. "Kinetic" momentum- 
energy tensor terms 

lim — • Udi(p <9V) + (d t Xa d\ a ) 

n-»3 n K 
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except directly dependent on temperature items, 

1 p 2 N lP N_zl y2 P2N2P 
3 p e lP 3 p e 2 p 

contain divergent vacuum terms, 



n \ p 2V n e lP p 2V n £ 2 p 

which should be renormalized. The parameter of dimensional regularization method is intro- 
duced via volume V n in non-integer dimension space, which is connected with integer-dimension 
space volume by simple relation: 

v n = Ni n = m 3 ■ r 2£ = vr 2e . 

Renormalization of divergent terms leads to the following result: 

Substitution of results of renormalization of divergent terms and two-point WF, expressed via 
effective masses and order parameter, to @ gives: 

= - l -[J 2 { T)mi {T)) + J 2 (T,m 2 (T))} + U( mi (T),m 2 (T),v(T)) (26) 



Free energy density ( |2T)| ) is distinguished into two parts: the first one describes contribution of 
ideal gas (temperature integrals over Bose-Einstein distribution), and the second one contri- 
bution of so called "self-consistent field" (indirectly dependent on temperature). 

It follows from fl24|) and (p6|), that generating functional, being considered on solution of gap 
equations for effective masses and equation of state for order parameter, turns into equilibrium 
free energy density. 

F(T, mi (T),m 2 (T),v(T)) = (27) 

2, ,4 



For this reason arbitrary constant in (f2~3| ) was chosen to be equal to 



N 2 fi 



16A(iV + 2)' 



Thermodynamic observables. Specific heat capacity, entropy and sonic speed are obtained 
from equilibrium free energy by usual thermodynamic relations. 



,2 



S 



V~~V\df) v ,Cv ~^V~ V\ <9T 2 ) V ,U c v 



Due to mentioned property of generating functional (|27|) and in view of (]23|) are identities on 
solution of gap equations and equation of state 

1 (dT{T,V)\ (dF\ T ( d 2 T{T, V) \ 

s ~~v\ df L~~[df ' Cy ~~v[ df° 1 

\ / V \ / v.mi ,mo. \ 



-T 



v,m\,m2 \ / V 

d 2 F\ ( d 2 F \ dv ( d 2 F \ dm x ( d 2 F \ dm : 

df 1 ) + [ dTdv ) df + [ dm x dT ) ~df + [ dTdm 2 j If 
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Effective masses and condensate order parameter derivatives are obtained by derivation of ( |21| ) 
and fl22j) , considered as identities. Result of calculation is given below. Specific heat capacity 
in low-symmetry phase: 



c v = - (4[J 2 (mi,T) + (iV-l)J 2 (m 2 ,T)] + 5 mp^, T) + {N - l)m\j x (m 2 , T) 



+ 



+ 



mlJ (mi, T) + (N - l)mjj (m 2 , T) 

4A(iV + 2 



-ti(m 2 ,T)(N-l] 



N 



4A(iV + 2) 
/ 1 (m 1 ,T) + l 



/ 2 2 (m!,T) 



' 4A(iV + 2) 

3 jj fi{m 2 ,T) 



+ 2{N-l)f 2 (m 1 ,T)f 2 {m 2 ,T) 



Here following designations are introduced: 

1 



fi(m,T) 



8tt 2 



In 



J (m, T), / 2 (m, T) = 2 Ji(m, T) + m J (m, T), 



, s * _ ! _ p^>/ l(m2 ,r) - 3 ) («&tS/ l(mi ,r) + 1). 



Specific heat capacity in high-symmetry phase: 



cv 



N 



4J 2 (m, T) + 5m 2 Ji(m, T) + m 4 J (m, T) — 



fl(m,T) 



V 



fi{m,T)- 



N 



2\{N + 2)/ 

Specific entropy in high-symmetry phase is easily obtained from that in low-symmetry phase: 

S = T it [ J2(mi ' T) + (iV - l )H m ^T)\ + [mlJ^muT) + (N - l)^^^)]) 

Phases stability conditions. Free energy density (p6|), obtained earlier by clearly thermo- 
dynamic relation (|25|) , can be calculated via statistical sum: 



-Tin Z. 



(28) 



The system energy spectrum, which is needed for the calculation, is determined by infill num- 
bers, effective masses and condensate values: 

E = Y. £ ipNip + (N - l)J2 £ 2k N 2k + V ■U(m 1 ,m 2 ,v), 



P 



k 



where Nip and N 2 p are numbers of if and \ a particles respectively. For the theory to be 
self-consistent in approximation used, one should sum over infill numbers Nip and N 2 p, but 
not the effective masses and order parameter, considering m 1; m 2 and v as constants (as it was 
stated in 0, §30]). After summation statistical sum is represented as a product of two terms: 
the first one corresponds to ideal gas and the second one to "self-consistent" field. Taking 
the logarithm and summing over impulses gives a function which coincide with generating 
functional. Consideration of gap equations for effective masses (0) and equation of state for 
order parameter ( pT] ) leads to the specific equilibrium free energy density (|26|). 

Keeping symmetrical parameters as constants, when calculating statistical sum, corresponds 
to the algorithm of calculation of non-equilibrium Landau statistical sum Zl M, which allows 
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to obtain non-equilibrium free energy functional Tl- Since the only symmetrical parameter in 
our case is condensate v, to obtain non-equilibrium Landau functional, one should substitute 
effective masses on order parameter dependence into generating functional. 

The functional obtained F(T,v) = F(T,mi(T,v),m2(T,v),v) reproduce equation of state 
for order parameter: 

d ( OF \ dm, f dF\ dm 2 fdF\ f dF\ 

\ 1 / T,m,2,v \ z / T,mi,v \ / T,mi,m,2 \ / T,m\,m,2 

and on solution of that equation v = v(T) it turns to free energy density, i.e. it possess all 
properties of non-equilibrium functional. 

A phase stability condition is order parameter Landau functional second derivation positiv- 

ity: 

d 2 F 9 24A , dmx , , 

— - = m\ - —v 2 + 2 m lV —±. 29 

dv z JM dv 

Inclusion of formal external classical source, interacting with condensate, in Lagrangian changes 
equation of state for 

C > C + p a (j) a , Pa = S a , N -V 

order parameter, which now depends on external source v = v(rj). In Hartree-Fock approxi- 
mation: 



2 n 

m! v 

1 N 



7]. 



Generalized sensitivity at rj — > turns out to be in inverse proportion with stability condition 



'dv\ (d 2 F\ 1 



\ dr lJ v ^o \ d v 2 J 

Being usual in Landau phase transitions theory, this result proves correctness of non- 
equilibrium free energy and stability condition fl29"|). 



Appearing in (29) condensate derivative of mass is obtained by order parameter differ- 
entiating of gap equations for effective masses (E3), considered as identities. In particular, 
low-symmetry phase stability condition looks like: 

(4A<* + 2) _ A fi^) /l(mi , r) + A 



dv 2 



2m{ 



N " K " ' \ N 



N-l \-N+l 



In the vicinity of critical temperature T C2 stability condition denominator is a slow variable 
function. Numerator temperature derivative, which designation was introduced earlier, 5 T ~ 

<5 _1 , this implies that 5 ~ \JT C2 ~T at the vicinity of T C2 . 

Specific heat capacity of low-symmetry phase contains S in inverse proportion, thus it is 
formally divergent at T C2 temperature. Sonic speed, which is inverse proportion with heat 
capacity, tends to zero. 

Stability condition of high-symmetry phase proportional to effective mass squared 

d 2 F _ 2 
dv 2 

which allows to obtain value of critical temperature T Cl analytically: 

_ / 3iV 
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Numerical calculations in Hartree— Fock approximation. Results of numerical calcu- 
lations of effective masses, order parameter and thermodynamical observables on temperature 
dependence for N = 2, iV = 4 and N = 9 and A = 0.01 ■ 2tt 2 are given below. 

If parameter A is chosen to be A = mi vac , mi vac = m\{T — 0), gap equations for effective 
masses and equation for order parameter posses a solution m\ = 2/i 2 , m| = 0, so that quanta 
of Xa fields turn to Goldstone bosons. 

It is convenient to normalize effective masses, order parameter and observables to m x vac . 

T 

Graphs of their temperature dependence are given in mentioned normalization (T =>- 

and so forth.). 




Fig.l Temperature dependence of effective masses, order parameter, sonic speed and specific 

heat capacity in Hartree-Fock approximation. 

Equations ([22|) and fl2T|) in low-symmetry phase admit thermodynamically stable in tem- 
perature region .. T C2 (solid lines at figures l.(a) and l.(b)) and in temperature region T C1 .. 
T C2 thermodynamically unstable (not shown) branches of solutions. In high-symmetry phase 
only one, stable, branch of solutions exists (dashed line). 

Below temperature T Cl in high-symmetry phase and above T C2 in low-symmetry phase, 
where stability conditions of the phases turn to zero, equations for effective masses (p2|) and 



order parameter (^1|) admit no solutions. When calculating correction to Hartree-Fock approx- 
imation, this criterion allows to easily find phases stability thresholds. 

As it follows from order parameter on temperature dependence, at phase equilibrium point 
T Cl < T c < T C2 , where equality of free energies of both phases is achieved, order parameter is 
small but not zero, consequently for finite N first type (close to second type) phase transition 
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takes place. Decreasing of condensate value at equilibrium point T c and decreasing of phase 

AT T c — T c 

coexistence region to T c ratio £ = = C2 — — with increase A/" (£|at=2 ~ 0.021, £|at=4 ~ 

0.011, £|tv=9 ~ 0.006) indicates that at the limit iV — > oo phase transition of second type takes 
place. Due to zero \a fields quanta masses, sonic speed is not zero at zero temperature. 

Numerical calculations confirm mentioned above conclusion of heat capacity divergence at 
vicinity of T C2 . However, it is well known, that heavy fluctuations in this region "blur" jump 
of heat capacity and other thermodynamic observables. 

Corrections to Hartree-Fock approximation turn out large at phases equilibrium point, 
which confirms inadaptability of the approximation at this temperature and, consequently, 
formality of made above conclusion about heat capacity divergence. 



5 Corrections to Hartree— Fock approximation. 

Calculation of corrections to Hartree-Fock approximation, which implies calculation of higher- 
rank WF, is important not only for estimation of the Hartree-Fock approximation temperature 
region of adaptability, but also for proof of correctness of mentioned in § renormalization 
method. 



Algorithm of calculations. Since equation for order parameter (0), which can be rewritten 
as follows 

2 8At> 2 4A / \ 

m i — + "T7 [ C 03(xxx) + (A/ - 1)C 21 ( XX \ X ) J 



0. 



(30) 



N N 

contains three-point WF, at first step one should derive partial solution of equation for WF 
C 03 ( XlX2X3 ) and C 21 concerned with sources. Proximate system of equations for three-point WF 
is derived from ([13]), (0) and ( ]T5|) by neglecting of higher-rank WF: 



^^ X G 03 ( X Xl X2 ) + m 1 C oz (x Xl X2 ) 



+ m\ C 2 



(ill \X2 ) 



24A 

'Iv" 

8A 
"iV 



G(x X1 )G(xx 2 ) + T^A<P(xxi )Atp 

(X X'2 ) 



D(x2X\)D(x2X) ~l~ -yy "X.{x2 x ) 3(( x 2 x l ) 

A 



^H^xG 2 l(xXl \x%) m 2 G 2 l(xXl |aJ2 ) ~ ^{x X\ )G( X X2 ) 



(31) 

(32) 

(33) 
(131 



In the framework of the iterative procedure used, for solving equations ([H]) and ( p2| 
two-point WF G ( XXl ) and D Q ^ XXl ) calculated in Hartree-Fock approximation are used. 

As accurate within first order of vanishing terms Xv 2 = ~~^ m i ? corrections to WF G^ XXl ) 
and D x t xxi \ are of the first order of vanishing and should be taken into account along with 
G 03 ( X1 X2X 3 ) and C 21 ( XlX2 \ X3 ). Thus, at second step equations for G^ XX1 ) and D^ XX1 ), which 
derived from (|TTD and (^) by substitution Xv 2 , expressed through mf, substitution of three- 
point WF and neglecting of the second order of vanishing sources (in particular four-point WF, 
multiplied by A). 



d^xG^xxi) + r m\ G^ XX1 ) 



7717 



3C 03 ( XXXl ) + (N — ljC^^a;!^) 



(34) 



d^D^^) + ml D l{xXl ) ~ -ml C 2l{xxi \ x) . (35) 

Taking into account tree-point WF and corrections to two-point WF in equation for order 
parameter and equations for effective masses leads to new temperature dependence of these 
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values in compare with Hartree-Fock approximation. New dependencies are found at third 
step. 

Preliminary conclusion of corrections to mi, m 2 and v magnitudes one could made with- 
out numerical analysis. Assuming the corrections small, introduce effective masses and order 
parameter as follows: 



rrii 



m\ + 5m 2 , m l = m lo + v 2 = Vq + 5v 2 



(36) 



Here mio, m 2 o and v corresponds to self consistent field approximation. Substitution of (|3q ) 
in (|30|) and @, (10) and separation by orders of vanishing gives the following expressions for 
corrections: 



5m\ 



4A(iV + 2) 



2(N-l)D i 



(xx) 



(^) /l(m20 , r) _3 



N 



+ (N- 1)C 2 

\{xx \x ) J 



(37) 



5m\ 



4\{N + 2) 



8 ■ N 



-2D, 



[xx) y 



fAX(N + 2) 



N 



A(mi ,T) + l +2G l(M) + 



Smf ■ N 



(38) 



^ = + 5 > + (iV " ^^(-l-)) " ( 39 ) 

Designations fi(m,T), which up to a coefficient is a squared mass derivative of two-point WF 
considered in Hartree-Fock approximation, 

dD o{xx) 



fi(m 10 ,T) 



(9m 2 
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T C2 — T, corrections, 



and 8 were introduced in chapter As long as at the vicinity of T C2 8 <• 
calculated by fl3"T|), ( |3"8"D and (|3TS|), are of high magnitude in this region. 

At fourth step corrections to thermodynamical observables are calculated. Neglecting of 
four-point Wightman functions, which contributions are of second order of vanishing, and 
substitution of two- and three-point WF, expressed through effective masses and condensate 
value, in "potential" terms of momentum-energy tensor gives: 



(T;> = hm (d, {x) dl Xi) G (xxi) + (N- l)d^ x) &( Xi) D {xxi 



N 



i\(N + 2) 



(N + l)m\ + 3(N - \)m\ N 



+/i 2 (m 2 + (N - l)m 2 ) + 



2 2 i 

m x v + 



2 

6Aw 4 



2 2 , 

m 1 m 2 + 



(40) 



N 



Since in high-symmetry phase three-point WF are identically equal to zero, corrections are of 
the second order of vanishing, and should not be considered in approximation which takes into 
account only effects of zero and first order of vanishing. 



Calculation of three— point WF. In Hartree-Fock approximation two-point WF G( XX1 ) 
and Dr xxi ), and commutators look like: 

-Sp Kp + 2 ( ip(x-xi) , ip(x 1 -x)\ n n k + 2 ( ik{x-xi) , ik(x 1 -x)\ 
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e ip(x-xi) _ e %p(xi-x) 



k 2£ k V 



gik(x —xi ) gifc(xi —x ) 



so that sources in equations ([H]), (|52"D, Q55D and Q33D, ( ff5j ) are sum of exponents. Since exponent 
is an eigen-function of differentiation operator, action of D'Alamber operator at each terms re- 
duce to multiplying by correspondent function of masses and impulses f(q, m). Hence, solutions 
of the equations are functions which differ from sources by factors f{q 1 m)~ l of exponents. 

Use of this algorithm and consideration of symmetrical properties of Wightman functions 
C 03 (xx lX2 ) gives: 



Re 



12X 2 n r d n p d n p' 



«->3 jV(2vr) 2n J e p e % 



cxp ■ 



12; Vexp^-1 ' 2 



e ip(x-xi) e «p' (x-x 2 ) _|_ e ip(xi-x) e ip' (x 1 -x 2 ) _|_ e %p{x 2 -x) gip' (x 2 -xi) 



m l - (P + P') 



(41) 



+ 
1 

"12 



e ip(x-xi) e ip'(x 2 -x) _|_ e ip(x!-x) e ip'(x 2 -xi) _|_ gip(22— x) gip' (X1—X2) ' 



m 2 — (P — P') 2 



+ 



Re 



' e ip(x-xi) e ip'{x-x 2 ) _|_ e ip(xi-x) e ip' (xi-x 2 ) _|_ e ip(x 2 -x) gip' (x 2 -x±) 

m i ~ (p + p') 2 



e ip(x—xi) e ip'(x 2 -x) _|_ gip(xi-x)gip'(x2-xi) _|_ gip(x2-a;)gip'(xi-X2) ' 



m 2 — (p — p') 2 

Here p = (p°,p) = (yp 2 + m 2 ,p), p' = (p' ,p') = (/p' 2 + m\ , p'). To obtain (||) G o(a;a;i) , 
D ( X xi)> Aip ( XXl ) and Xl ) one should turn from summation to integration over non-integer 

dimension space. 

In solutions of (|32|) and (|33|) terms — ^ —7 and — 5 —77 respectively arise. 

wij — (p — k) 2 mf — (k + k') 2 

Denominators of these terms turn to zero for some impulses p, k if condition mi > 2m 2 is 
satisfied. Physically satisfaction of this condition corresponds to possibility of ip to pair Xa 
decay. In tree approximation decay width is: 



N - 1 2A V 



N 2 irnii 



1- 



2m 2 
mi 



Let us now introduce helper function C^ xxi u a %, which satisfy to equations obtained from (]32|) 
and (|33| ) respectively by substitutions 

+m 2 i^ d^ 2 + mj- imiT, 0„fi£ + m 2 , =^ <9 M d£ + m 2 - im 2 R 



Since WF are real, by definition C 21 ( XXl \ X2 ) = - 



au _j_ au * 

2i(a:a:i |*2 ) ' i\{xx\ \x 2 ) 



After consideration of C 21 ( XXl \ X2 ) symmetries with regard to arguments permutation: 



^C 21 ( XXl \ X2 ) — lim 



-4A? 



3 N(2ir] 



in 



Re 



'Re ( { 



d n p d n k 
£p£k V^exp^-1 2 



exp — 1 2 



x 



ip(x— xi)gik(x— X2) _|_ gip(xi— x)gik(xi— X2) 



+ 



m i — (p + k) 2 — im 2 T 

J^g ^gjp(x-xi)gik(x2-x') _|_ gip(xi-x)gik(z2-i'i) 



+ 



m 2 — (p — k) 2 — zm 2 r 



(42) 
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— lim 



-4A 



N(2ll 



+ 



12 



Re 



Re 



Lf 

• 6 J 



Re ^■H^-x) e ik' (x2-xi)^j ^ e ik(x-X2) e ik' (X2-X1) 



d n k d n k' 




m 2 — (k + k') 2 — im\Y 

Re ( e ik ( x 2-x) e ik' (x 2 -x!) 



(k-k 



A2 



J^g / e ik(x-x 2 ) e ik' (x 2 ~xi) 



m 2 — (k + k') 2 — imiT m 2 — (k — k') 2 — wiiT 




Here k = (A; , k) = (^Jk 2 + to?, k), k' = (A;' , #) = (^k' 2 + to 2 , k'). Solutions, written above, 
are symmetrical of permutations of arguments. Taken in coincide points, they depend on inner 
system parameters, but not 4-coordinates. 



Wightman functions C 03 f xxx ) and C, 



2i(xx\x) automatically separate into items of three types. 
Items of the first type contain product of Bose-Einstein distributions under the integral and 
are finite at the limit n — > 3. Items of the second type ("mixed integrals") which contain 
under integral Bose-Einstein distribution for one integration variable and items of the third 
type ("vacuum integrals") which don't contain Bose-Einstein distribution under integral are 
divergent at the limit n — > 3. 

Asymptotic expansion of integrands is used for regularization of divergent terms of three- 
point WF in the framework of dimensional regularization. Three-dimensional impulse p (k) to 
effective mass TOi (to 2 ) quotient serves as expansion parameter, so that region of integration is 
divided into four regions where values of these parameters are greater (less) than one. In each 
region an expansion valid in the region is used. 

The integrals are well defined in space of dimension 1 < n < 2. Since "vacuum integrals" 
are divergent by both integration variables, analytical continuation into n = 3 dimension space 

lead not only to e^ 1 pole (which is responsible for In I 1 terms in renormalized integrals), 



but also e 2 pole (which is responsible for In 2 




terms) . 



As "mixed integrals" are divergent only by one integration variable, at the limit n — > 3 only 



/ 771 \ 

e^ 1 pole arise. Coefficient at In I \ , formed by Bose-Einstein distribution, is a function of 

particle mass and temperature. 

After poles subtraction one obtain the following expressions for three-point Wightman 
functions: 

-18A 



03(iC X X )ren 



21 (:£ X \ x )ren 

mi m 2 



Q 

-21n, v ., 

A 2 

+C 3 m 2 In 



/i°(m 1)mi) mi , T) + J^rm, T) (c t - 2 In ^ + 

+K° 1 (m 1 ,m 1 ,m 1 ,T) + C 2 m 2 In 2 (^j + C 3 m 2 In (^j + C 4 m 2 , 

/ 1 °(mi,m2,m 2 ,r) + i ^(to^T) + ^|Ji(m 2 ,T)^ (Ci- 



(2ir) 2 N 



-4A 



(2tt) 2 N 

+ ^K°(m 1 ,m 2 ,rn 2 ,T) + ^K°(m 2 ,m 1 ,m 2 ,T) + C 2 m\ In 2 ' ' l '" 2 



77li7n 2 



+ C 5 m 2 + L°(toi,to 2 ,to 2 , 1) 



+ 



-2A 



(2tt) 2 N 



+K 1 °(m 2 ,m 2 ,m 1 ,T)+ (2-^|j Jx(m 2 ,T) fd 




V A 2 
I^(m 2 ,m 2 ,m 1 ,T)+ 

- C 5 m 2 + C 6 m 2 + 



+ 
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+C 2 (2m 2 ~ m\) In 2 (§) + C 3 (2m 2 - m\) In ($) + L°(m 2 , m 2 , rm, 4/3) 

Coefficients Cj were calculated approximately: Ci « —1.386, C 2 ~ 0.044, C 3 0.013, C 4 « 
-0.031, C 5 « -0.012, C 5 « -0.031, C 6 « -0.089. Introduced above functions I n (x,y,a,T), 
K n (x, y, a, T) and L n (x, y, a, /3, T) are defined as follows: 

2 r rd 3 pd 3 k 1 1 



•Re 



(2vr) 4 
(p + k) 2m 



e p e k 
+ - 



£ P i £ k 
exp - 1 exp 



(p-k) 



2m 



a 2 — (p + k) 2 — iaT) n (a 2 — (p — k) 2 — iaT) r 



(43) 



K^{x,y,a,T) 
•Re 



|p|=oo r \k\=vd 3 pd 3 k 



(2?r)- j pi i. 
p+ fc) 2m 



|Jfe|=o £p£fe exp ^_l 
(p - fc) 2m 



a 2 - (p + k) 2 - zoT) n + (a 2 - (p - k) 2 - zaT) r 
1 /-[Pl=* f\ k \=yd 3 pd 3 k 



(44) 



•Re 



L™(x,y,a,/3) = 

2(27r) 4 7|p|=o 7 |*|=o £p£fc 
/?(p+fc) 2m (2-/5)(p-*) 2m 



+ 



(45) 



(a 2 — (p + k) 2 — iaT) n (a 2 — (p — k) 2 — iaT) 
In fl43D , (|44]) and d45|) 4-impulses p and k depend on arguments x and y as: p = (p°,p) = 

(ep, p) = (vW^ 2 , p), k = (*°, *) = (e k , k) = (v/fc^T^, *). 

Proportional to introduced in chapter |] temperature integral Ji(m,T) terms are renor- 
malized "mixed integrals". Renormalized "vacuum integrals" are not explicitly temperature 
dependent. 

( jjl 2 \ / ^ 777,2 \ 

— | j and m\ In ( — — — J is 

logarithmically divergent at the limit T — > . This divergence, caused by long-range action, is 
described by diagrams of massless Xa quanta exchange. Since at the most interesting region of 
phase transition m 2 is not zero, this problem was not investigated in details. 



Calculation of corrections to two— point Wightman functions. Taking into account 
fliHD , (H|]) and fljjD it is convenient to rewrite equations fl34]) and ( |35"D as follows: 



+ m 2 } [d^ + ml 
AXmj(N 
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36A m\ r 1 
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(46) 
(47) 



Since, like in the case of three-point WF, divergent for mi > 2m 2 and some impulses terms 
1 1 

and - — n - ; arise in solutions, equations for two-point WF should 



(ml - (p - k) 



2^2 



(ml - (k + k') 2 ) 2 



be redefined. Use of the same way of redefinition as was used for C 21 



(x x\ |xa ) 



leads to: 



(xx 1 ) 
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1 V 
exp^-l + 2, 
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Refe*^ 1 -x) e ik'(xi -x)\ He(e ik ( Xl ~ x ) e ik '( x - x ^)^ 
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AA(in) = lim 
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n^3iV(27r) 2 ™y Uxp^-1 2 
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Re 



j^gjg^p^i -x) e ik(xi -x)^ Re(e ip ^ Xl - x ) e i]i ( x ~ x i )) 



(49) 



+ 



(m 2 — (p + k) 2 — im 2 T) 2 (m\ — (p — k) 2 — im 2 T) 2 ) 

After renormalization of corrections to two-point WF, considered in coincide points, in the 
framework of dimensional regularization one obtain: 



G 1 (xx) r 



m\ ll(m u m^m x ,T) + h(m x ,T) ^-21n^^ + 
+m\ K° 2 ( mi ,m 1} m 1} T) + C 2 m\ In 2 + C 7 m\ In (^j + C s mj 
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Ji(mi,T) H ^J 2 (m 2 ,T) 
m 2 J 
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YK 2 (rn 1 ,m 2 ,m 2 ,T) + -^Kl(m 2 ,mi,m 2 ,T) + C 2 m 2 ln 



CV _ 
2 

2 2 (m 1 m 2 



A 2 
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+C 7 m 2 In 



mi m 2 



Cn m\ + L 2 (mi, m 2 , m 2 , 1) 



Here C 7 « -0.568, C 8 w -0.929, C 9 » -0.111, Ci » -0.818, Cn w -1.042. 

In order to calculate (specific) equilibrium free energy, which is expressed via momentum- 
energy tensor components by relation (|25|), considering correction to Hartree-Fock approxima- 
tion of the first order of vanishing, one should derive and renormalize derivatives of corrections 
to two-point WF. 

Differentiation of (|5]) and (^) gives respectively: 

Vn ' di ^ d t x i) G ^ xx ^ = 
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After renormalization one obtain: 
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Numerical values of the coefficients C, are: C 12 ~ 0.719, Ci 3 = 0.011, Ci 4 ~ 0.177, Ci 5 ~ 
-0.019, Cie « 0.065, CiT « -0.054, Ci 8 « -0.021, Ci 9 « 0.198, C 20 « 0.026, C 21 « -0.046, 
C 22 « -0.526, C 23 « 0.617, C 24 « 0.168, C 25 « 1.341, C 26 « 0.396, C 27 « -0.636. 
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Thus specific free energy density calculated accurate to first order of vanishing terms looks 
like: 



mlv 



4Xv 4 



lim G" (v r , +(N-1) lim £>", TT H - 

where ( |2~3| ) is evaluated on solution of flTU|) and (^0). 



+ F(T, mi (T),m 2 (T),v(T)), 



Numerical calculations and conclusions. Taking into account correction to Hartree-Fock 
approximation does not bring to appearance of a new branch of solutions of equations for ef- 
fective masses and condensate. These values on temperature dependencies are depicted at Fig. 
2 (solid lines at Fig 2. a, Fig 2.b, Fig 2.c). For easier comparison, by dashed line these values 
on temperature dependencies in Hartree-Fock approximation are depicted. 




Fig. 2 Effective masses, order parameter, equilibrium free energy and specific heat capacity on 
temperature dependencies calculated accurate to corrections of the first order of vanishing. 

From Fig. 2 it follows that in low-symmetry phase only one branch of solutions remain 
thermodynamically stable in temperature region < T < T C2 . 

Significant deviation of m<i on temperature dependence in compare with that in Hartree- 

(771^ \ 
) 

and m 1 in — -= — 

V J\ 



in three-point WF and corrections to two-point WF. 
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Critical temperature T C2 decreases in compare with that calculated in Hartree-Fock ap- 
proximation. With increase of N, relative magnitude of corrections, estimated as quotient 

7 = Tc2 P y - Tc2 ■ 100% decreases (7^=2 « 3.26%, 7^=4 « 3.21%, 7^=9 ~ 2.61%) and at 

the limit N ^ 00 tends to zero, which agrees with conclusion made in chapter 

Though 1712 on temperature dependence significantly changes, at low temperatures thermo- 
dynamical observables (equilibrium free energy and heat capacity) values are close to those 
in Hartree-Fock approximation, as expected from general conclusions. The closer to T C2 , the 
greater deviation of observables from that in Hartree-Fock approximation, in particular heat 
capacity, which agrees with made above conclusion of formality of heat capacity divergence at 
T C2 ( ) vicinity. 
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